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National Phase of: PCT/CH99/00476 
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Assistant Commissioner for Patents 
Washington, DC 20231 

Sir: 

PRFT.TMTNARY AME NDMENT 

Prior to calculation of the filing fee and in order to place the above identified application 
in better condition for examination, please amend the claims as follows: 
IN THE SPECIFICATION 

On page 1, after the title, please insert --This application is a 371 of PCT/CH99/00476.- 

TN THE CLAIMS (AS AMENDED) 

Please substitute the following amended claims for corresponding claims previously 
presented. A copy of the amended claims showing current revisions is attached. 

1. (Amended) A method for local and superficial characterization of a turbid 
medium using the following parameters: 

1) the refractive index n of the turbid medium 

2) the absorption coefficient fi^ of the turbid medium 

3) the reduced scattering coefficient pi/ of the turbid medium 
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4) the phase function parameter Y-{l-g^)/{l-g^) of the turbid media, where g^and are the first 
two moments of the Legendre polynomial development of the phase function p, 
and comprising the steps of: 

- measuring the spatially-resolved reflectance R(p }of the turbid medium (p being the source- 
detector distance) by any means, comprising an illumination beam as a source and an optical 
detector, which, by using optional signal processing, which may involve filtering and de- 
convolution operation to correct for the non-zero area of either the illumination source or the 
detector, allows for the precise determination of the said spatially-resolved reflectance R{p), 

- mathematically processing R{p) to compute at least one of the said parameters: n, /i^ , jLt/, 7 
and/or the variations, in time and/or space, of at least one of the said parameters: hn, /S^fi^ A jit/. 
Ay, whereby an "inverse problem", which consists in extracting the optical coefficients from the 
spatially resolved reflectance data is solved, and whereby a "direct problem" consists in 
computing the spatially resolved reflectance from the values of the optical coefficients n, ji^ , /x/, 
7 involved in a model of propagation of the light in turbid medium and whereby the said "model" 
incorporates a Legendre polynomial development to the second order of the said "phase 
function", and whereby the said "phase parameter" y is introduced in the computation as an 
independent parameter. 

2. The method of claim 1, wherein said spatially resolved reflectance is measured by 
a probe comprising at least one optical fiber carrying the light from the source unit to the turbid 
medium and at least one optical fiber collecting the reflected light and carrying it to detection 
unit, whereby the combination of a variety of emitting fibers and of receiving fibers yields a set 
of distances p at which the reflectance R(p) is measured. 
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3. (Amended) The method according to claim 1, wherein said spatially resolved 
reflectance R(p) is measured for a set of values of p, by using a probe composed of optical fibers 
in any of the following configurations: 

one emitting optical fiber and a set of optical receiving fibers 

a set of optical emitting fibers and one optical receiving fiber 

a set of optical emitting fibers and a set of optical receiving fibers 
giving the spatially resolved reflectance R(p) at a variety of source-detector distances p and 
wherein the emitting and receiving fibers are arranged in one of the following configurations: 

on a line, 

on crossed lines, 

on a circle, 

on an ellipse 

on crossed ellipses 

on a disk, a rectangle, or any surface, as a dense and contiguous arrangement of fibers, 
on any pattern resulting from the combination of the above mentioned patterns. 

4. (Amended) The method according to claim 3, wherein said spatially resolved 
reflectance is measured by an optical and electronic micro-system comprising a collimated or 
focused beam as illuminating source and ID or 2D arrays of optical detectors. 

5. (Amended) The method according to claim 4, wherein either the probe or the 
optical and electronic micro-system are put in contact to the turbid medium. 
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6. (Amended) The method of claim 1, wherein said spatially resolved reflectance 
is measured by a non-contact system, comprising at least one of the following combination of 
optical systems: 

a fixed optical system to irradiate the turbid medium with a collimated or focused beam 
forming the illuminating source and a fixed optical system comprising an imaging system 
forming the image of the measured area of the turbid medium on a said "optical detector", which 

2 can be formed of ID or 2D array of optical detectors, whereby this second optical system can be 
5 identical to the first one and whereby the array of optical detectors can be either one of the 

following systems: 
■ a set of optical fibers, 
p ■ an optical and electronic micro-system (MOEM), 

b - a fixed optical system for the collimated beam illuminating source and a scanned optical 

3 system for the said "optical detector", 

a scanning optical system for the collimated beam illuminating source and a fixed optical 
system for the said "optical detector", 

a scanning optical system for the collimated or focused beam used as an illuminating 
source and a scanning optical system for the said "optical detector". 

7. (Amended) The method according to claim 1 , wherein the absorption 
coefficient fi^ the reduced scattering coefficient ft/ and the phase function parameter Y are 
determined by fitting the measured spatially-resolved reflectance R(p, /J,/, fi^ j) to a set of 
discretized data obtained by using Monte Cario simulations, or to interpolating functions giving 
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continuous approximation of the discretized data obtained by Monte Carlo simulations, and 
whereby said "Monte Carlo simulations" are based on a photon propagation model comprising a 
phase function approximated by a Legendre polynomial development Hmited to the second 
order. 

8. (Amended) The method according to claim 7, wherein one or more of the 
following signal processing steps are performed: 

fitting the measured reflectance R(p) hy the function: 
m^p'^ exp(m3p) 

to give the values of the parameters m^ m^and m^, assuming that the expression 
R(p)=myp'^ exp(m3p) gives a smoothed expression of the spatially resolved reflectance R(p), 

computing the slopes — and — (in (p)) , or any mathematical combinations 

of these two latter quantities and R(p), from analytical functions using the parameters mp m^ m3, 
or by numerical procedures from the expression R(p)== m^p"^ exp(/?^p) , 

computing the values of at least one of the said parameters: n, ix^ , /x/, y and/or the 
variations, in time and/or space, of at least one of the said parameters: An, A/a^, A /x/, Ay from the 
relationship between the reflectance intensity R(p) and the slope of \nR(p) (denoted d p InR), 
determined at a fixed distance p comparable to the transport mean free path, whereby the 
computation is made from the data obtained by Monte Carlo simulations, provided that said 
"Monte Carlo simulations" are based on a photon propagation model comprising a phase 
function approximated by a Legendre polynomial development limited to the second order. 
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9. (Amended) The method according to claim 7, wherein the following signal 
processing steps are performed: 

computing the reduced scattering coefficient ' and the phase function parameter y by 
using the following form of the reflectance; 

R(p) = (A(p,r,ii:) +B(ix^,ii;)f 

where the function A(p, y, ft'j and B(n^, fij depend also on the sources and detectors 
characteristics, and the refractive index of the sample, and comprising the steps of: 
computing the slopes of the square root of the spatially resolved reflectance 

— jR{p.^ii,f.i^,Y) = — ( i"; > r ) . which depends weakly on the absorption coefficient fi^ for 
dp " " c)p 

0.3<p/x/<5, for at least two distances p, 

determining the parameters jLi/ and y by a polynomial interpolation of the data obtained 

by Monte Carlo simulations, whereby said "Monte Carlo simulations" are based on a photon 

propagation model comprising a phase function approximated by a Legendre polynomial 

development limited to the second order. 

10. (Amended) The method of claim 9, wherein the absorption coefficient is 
determined by using the equation: 



^/?(p,^;,/i,,y)-/(y,At;) 
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where / and h are continuous functions of the parameters p, fi/ and ythat can be approximated 
by a polynomial interpolation of the data obtained by Monte Carlo simulations, whereby said 
"Monte Carlo simulations" are based on a photon propagation model comprising a phase 
function approximated by a Legendre polynomial development limited to the second order. 

11. (Amended) The method according to claim 1, wherein the difference Afi^ =[1^- 
fi^ between the absorption coefficient ji^ and a known value fi^js computed from the quantity 

fil , ;W^o' r) ~yl^(P^ 1^1^ /^a^ Y) - B(a')- B(aJ), whereby the function B(a) of the albedo a' 
can be can be approximated by a polynomial interpolation of the data obtained by Monte Carlo 
simulations, whereby said "Monte Carlo simulations" are based on a photon propagation model 
comprising a phase function approximated by a Legendre polynomial development limited to the 
second order, and whereby the calculation can be done for a single y value, because the 
influence of the phase function parameter and /in B(a) are particularly weak. 

12. (Amended) The method according to claim 1, wherein the illuminating source 
is a broadband source or a white light source and the detector unit comprises a spectrograph or 
any wavelength analysis system to measure the wavelength dependence of at least one of the 
parameters (n, li^fi/, y) . 

13. (Amended) The method according to claim 1, wherein said turbid medium is a 
biological medium. 
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14. (Amended) The method according to claim 1 , wherein the measurement and 
processing is repeated at different locations of the sample, to build images of any one of the said 
parameters (n, y). 

15. (Amended) An apparatus using the method of claim 1 for local and superficial 
characterization of a turbid medium, 

a) comprising a source, a detection unit, reference means, signal processing means, a probe 
'% comprising at least one optical fiber connecting said source unit to the turbid medium and at least 
il one optical fiber connecting the turbid medium to the said detection unit, and reference means 
!B b) where the distance between the source and the detector is close to one transport mean free 

path. 

I'S 16. (Amended) An apparatus using the method of claim 1 for local and superficial 

characterization of a turbid medium, 

a) comprising an optical and electronic micro-system comprising at least one illuminating 
source, at least one detector, signal processing means and reference means, 

b) where the distance between the source and the detector is close to one transport mean free 
path 

17. (Amended) An apparatus using the method of claim 1 for local and superficial 
characterization of a turbid medium, 

a) comprising a coUimated or focused beam used as an illuminating source, at least an optical 
detector for the detection unit, a fixed or scanning optical system for the illuminating source and 
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a fixed or scanning optical system for the said "optical detector", signal processing means and 
reference means, 

b) where the distance between the source and the detector is close to one transport mean free 
path 

18. (Amended) A system comprising three apparatus described in claim 15, and 
characterized by the fact that the distance between the collimated or focused optical beam used 
as illuminating source and the emitting point connected to an optical detector varies from 0.1 to 
2mm. for application to biological media and to turbid media having a transport mean free path 
similar to biological media. 

19. (Amended) A test, where the control of the homogeneity of the sample over 
the probed area is performed with the apparatus of 15, which can be carried out according to the 
following procedure: in the apparatus, disposing two illuminating fibers symmetrically in regard 
to the collecting fibers; comparing the reflectance curves for each illuminating fiber to detect the 
heterogeneity of the investigated region or obstructions beneath the fibers; and, if the two curves 
are sufficiently close, validating the measurement and calculating the average of the two curves. 

20. (Amended) A calibration and normalization procedure, which are carried out 
with the apparatus of claim 15, whereby the following steps are performed: 

1) in order to perform relative intensity measurements, the differences of transmitted intensity 
between each fiber for the apparatus of claim 15 are determined by performing a measurement 
on a turbid phantom illuminated unifonnly or a diffusing sphere of perfectly uniform properties; 
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2) in order to perform absolute intensity measurements, a calibration performed on a turbid 
medium of known optical properties, which can be realized according to any one of the 
following recipes: 

a) solid or liquid turbid medium which properties have been measured by other techniques, 
or reported in the literature, 

b) water suspension of micro-spheres of known size distribution and refractive index. 



REMARKS 

Attached hereto is a marked-up version of the changes made to the claims by the current • 
amendment. The attached page is captioned " Version with markings to show changes made ." 
The above amendments are made to place the claims in a more traditional format. 

Respectfully submitted, 



NIXON & VANDERHYE P.C. 




J^hael J. Shea 
Reg. No. 34,725 



MJS:lmy 

1 100 North Glebe Road, 8th Floor 
Arlington, VA 22201-4714 
Telephone: (703) 816-4000 
Facsimile: (703) 816-4100 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 

1. (Amended) A method for local and superficial [(probed volume and measuring 
area with dimensions on the order of one transport mean free path)] characterization of a turbid 
medium using the following parameters: 

1) the refractive index n of the turbid medium 

2) the absorption coefficient [i^of the turbid medium 

3) the reduced scattering coefficient /x/ of the turbid medium 

4) the phase function parameter y =(l-g,)/(l-g,) of the turbid media, where g, and g, are the first 
two moments of the Legendre polynomial development of the phase function p, 

and comprising the steps of: 

- measuring the spatially-resolved reflectance R(p )of the turbid medium (p being the source- 
detector distance) by any means [mean], comprising an illumination beam as a source and an 
optical detector, which, by using optional signal processing, which may involve filtering and de- 
convolution operation to correct for the non-zero area of either the illumination source or the 
detector, allows for the precise determination of the said spatially-resolved reflectance R(p), 

- mathematically processing R(p) to compute at least one of the said parameters: n, , ii/, y 
and/or the variations, in time and/or space, of at least one of the said parameters: An, A^x^, A fi/, 
Ay, whereby an [the said] "inverse problem", which consists in extracting the optical 
coefficients from the spatially resolved reflectance data is solved, and whereby a [the said] 
"direct problem" consists in computing the spatially resolved reflectance from the values of the 
optical coefficients n, [i^ , n/, y involved in a model of propagation of the light in turbid medium 
and whereby the said "model" incorporates a Legendre polynomial development to the second 
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order of the said "phase function", and whereby the said "phase parameter" y is introduced in 
the computation as an independent parameter. 

2. The method of claim 1 , wherein said spatially resolved reflectance is measured by 
a probe comprising at least one optical fiber carrying the light from the source unit to the turbid 
medium and at least one optical fiber collecting the reflected light and carrying it to detection 
unit, whereby the combination of a variety of emitting fibers and of receiving fibers yields a set 
of distances p at which the reflectance R(p) is measured. 

3. (Amended) The method according to claim [any one of claims] 1 [and 2], 
wherein said spatially resolved reflectance R(p) is measured for a set of values of p, by using a 
probe composed of optical fibers in any of the following configurations: 

one emitting optical fiber and a set of optical receiving fibers 

a set of optical emitting fibers and one optical receiving fiber 

a set of optical emitting fibers and a set of optical receiving fibers 
giving the spatially resolved reflectance R(p) at a variety of source-detector distances [distance] 
p and wherein the [. The] emitting and receiving fibers are [can be] arranged in one of the 
following [a variety of] configurations [. In some particular arrangement they can be positioned]: 

on a line, 

on crossed lines, 

on a circle, 

on an ellipse 

on crossed ellipses 
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on a disk, a rectangle, or any surface, as a dense and contiguous arrangement of fibers[, 
such as it can be obtained with a bundle of fibers or a multi-core fiber], 

on any pattern resulting from the combination of the above mentioned patterns. 

4. (Amended) The method according to claim [any one of claims 1 , 2 and] 3 , 
wherein said spatially resolved reflectance is measured by an optical and electronic micro-system 
comprising a collimated or focused beam as illuminating source and ID or 2D arrays of optical 
detectors[, arranged in a way similar to the one described in claim 3, and which, in some 
embodiment, can be a CCD or MOS camera]. 

5. (Amended) The method according to claim [any one of claims 1 -]4, wherein 
either [where] the probe [of claims 2 and 3] or the optical and electronic micro-system [of claim 
4] are put in contact to the turbid medium. 

6. (Amended) The method of claim 1 , wherein said spatially resolved reflectance 
is measured by a non-contact system, comprising at least one of the following combination of 
optical systems: 

a fixed optical system to irradiate the mrbid medium with a collimated or focused beam 
forming the illuminating source and a fixed optical system comprising an imaging system 
forming the image of the measured area of the turbid medium on a said "optical detector", which 
can be formed of ID or 2D array of optical detectors, whereby this second optical system can be 
identical to the first one and whereby the array of optical detectors can be either one of the 
following systems [system]: 
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■ a set of optical fibers [arranged in a variety of configurations similar to those described 
in claim 3], 

■ an optical and electronic micro-system (MOEM)[, which in some embodiment can be a 
CCD or MOS camera, with a functionality similar to the one described in claim 4], 

a fixed optical system for the coUimated beam illuminating source and a scanned optical 
system for the said "optical detector", 

a scanning optical system for the coUimated beam illuminating source and a fixed optical 

system for the said "optical detector", 

a scanning optical system for the coUimated or focused beam used as an iUuminating 
source and a scanning optical system for the said "optical detector". 

7. (Amended) The method according to claim [any one of claims] 1 [. to 6.], 
wherein the absorption coefficient {i^ihe reduced scattering coefficient ix; and the phase function 
parameter y are determined by fitting the measured spatially-resolved reflectance R(p, ix/, /x^ y) 
to a set of discretized data obtained by using Monte Carlo simulations, or to interpolating 
functions giving a continuous approximation of the discretized data obtained by Monte Carlo 
simulations, and whereby said "Monte Carlo simulations" are based on a photon propagation 
model comprising a phase function approximated by a Legendre polynomial development 
limited to the second order[ In a preferred embodiment, cubic splines are used as interpolating 
functionsJi 

8. (Amended) The method according to claim [any one of claims 1 to] 7, wherein 
one or more of the following signal processing steps are [can be] performed [in addition to the 
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processing steps of any of the claims 1-6 whereby the processing of the spatially resolved 
reflectance data of claim 7 can be simplified and accelerated by]: 
fitting the measured reflectance R( p) by the function: 

m^p"^ txp[m^p) 

to give [This fit gives] the values of the parameters m„ m^ and m, , assuming that the expression 
R(p)= m^p"^ exp(m3p) gives a smoothed expression of the spatially resolved reflectance R(p)[.]^ 

i3 - computing [Computing] the slopes — ^/?(p) and — (in /? (p )) , or any mathematical 

0 combinations of these two latter quantities and R(p), from analytical functions using the 

Ui parameters m^ m^ m., or by numerical procedures from the expression R(p)= m^p"^ exp{m^p)^ 

- computing the values of at least one of the said parameters: n, [x^ , /a/, y and/or the 
JS variations, in time and/or space, of at least one of the said parameters: An, A/x,, A^/, Ay from 
13 the relationship between the reflectance intensity R{p) and the slope of InRfpj (denoted 

a p InR), determined at a fixed distance p comparable to the transport mean free path, whereby 
the computation is made from the data obtained by Monte Carlo simulations, provided that said 
"Monte Carlo simulations" are based on a photon propagation model comprising a phase 
function approximated by a Legendre polynomial development limited to the second order. 

9. (Amended) The method according to claim? [any one of claims 1 to 7], 
wherein the following signal processing steps are [can be] performed [in addition to the 
processing steps of any claim 1-6 whereby the processing of the spatially resolved reflectance 
data of claim 7 can be simplified and accelerated by]: 
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computing the reduced scattering coefficient ' and the phase function parameter y by 
using the following form of the reflectance: 

R(p) = (A(p.Y.li;) +B(li,li:)f 

where the function A(p, y, iij and B(n^, nj depend also on the sources and detectors 
characteristics, and the refractive index of the sample, and comprising the steps of: 

computing [Computing] the slopes of the square root of the spatially resolved reflectance 



— Ji?(p.;u;,/i„,y) = — ( p 4< ; , y ) , which depends weakly on the absorption coefficient fi^ for 

dp ' 9p 

0.3<p/y./<5, for at least two distances p, 

determining [Determining] the parameters /i/ and 7 by a polynomial interpolation of the 
data obtained by Monte Cario simulations, whereby said "Monte Carlo simulations" are based on 
a photon propagation model comprising a phase function approximated by a Legendre 
polynomial development limited to the second order. 

10. (Amended) The method of claim [claims] 9, wherein the absorption coefficient 
jLi^is determined by using the equation: 



^R{p,ixi,^„y)-f(y,fi;) 



where/and h are continuous functions of the parameters p, ]ij and y [of] that can be 
approximated by a polynomial interpolation of the data obtained by Monte Cario simulations, 
whereby said "Monte Carlo simulations" are based on a photon propagation model comprising a 
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phase function approximated by aLegendre polynomial development limited to the second 
order. 

11. (Amended) The method according to claim [any one of claims] 1 [to 7], 
wherein the difference Ajx^ =^'1^00 between the absorption coefficient /i^ and a known value iij^ 
computed from the quantity ^R(p,fi'„^,„Y) - ^R(p,fil,l^,,Y) = B(a ')- B(a;), whereby the 
function B(a') of the albedo a' can be can be approximated by a polynomial interpolation of the 
data obtained by Monte Carlo simulations, whereby said "Monte Carlo simulations" are based on 
a photon propagation model comprising a phase function approximated by a Legendre 
polynomial development limited to the second order, and whereby the calculation can be done 
for a single y value, because the influence of the phase function parameter and yin B(a') are 
particularly weak. 

12. (Amended) The method according to claim [any one of claims] 1 [to 1 1], 
wherein the illuminating source is a broadband source or a white light source and the detector 
unit comprises a spectrograph or any wavelength analysis system to measure the wavelength 
dependence of at least one of the parameters (n, ju.^ ju./, y) ■ 

13. (Amended) The method according to claim [any one of claims] 1 [to 12], 
wherein said turbid medium is a biological medium. 

14. (Amended) The method according to claim [any one of claims] 1 [to 13], 
wherein the measurement and processing is repeated at different locations [location] of the 
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sample [(multi-site measurements)], to build images of any one of the said parameters 
(n, Y)- 

15. (Amended) An apparatus using the method of claim 1 [to 3, 5 or 7 to 1 4] 
for local and superficial characterization of a turbid medium, 

a) comprising a source, a detection unit, reference means, signal processing means, a probe 
comprising at least one optical fiber connecting said source unit to the turbid medium and at least 
one optical fiber connecting the turbid medium to the said detection unit, and reference means 

b) where the distance between the source and the detector is close to one transport mean free 
pathi 

16. (Amended) An apparatus using the method of claim 1 [, 4, 5, or 7 to 14] 
for local and superficial characterization of a turbid medium, 

a) comprising an optical and electronic micro-system comprising at least one illuminating 
source, at least one detector, signal processing means and reference means, 

b) where the distance between the source and the detector is close to one transport 
mean free path 

17. (Amended) An apparatus using the method of claim 1 [or 6 to 1 4] for 
local and superficial characterization of a turbid medium, 

a) comprising a coUimated or focused beam used as an illuminating source, at least an optical 
detector for the detection unit, a fixed or scanning optical system for the illuminating source and 
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a fixed or scanning optical system for the said "optical detector", signal processing means and 
reference means, 

b) where the distance between the source and the detector is close to one transport mean free 
path 

18. (Amended) A system comprising three [Three] apparatus [similar to the 
apparatus] described in claim [any one of the claims] 15[, 16 and 17 ], and characterized by the 
fact that the distance between the collimated or focused optical beam used as illuminating source 
and the emitting point connected to an optical detector varies from 0.1 to 2mm. for appMcation to 
biological media and to turbid media having a transport mean free path similar to biological 
media. 

19. (Amended) A test, where the control of the homogeneity of the sample over 
the probed area is performed with the apparatus of [any one of the claims] 15[, 16,17 and 18 ], 
which can be carried out according to the following [optional] procedure: in [In] the apparatus^ 
disposing two illuminating fibers [are disposed] symmetrically in regard to the collecting fibersi 
comparing r. If the sample is homogeneous,] the reflectance curves for each [curve should be 
identical with either] illuminating fiber to detect the [. Therefore,] heterogeneity of the 
investigated region or obstructions beneath the fibers : and, if [are detected by comparing the two 
curves. If] the two curves are sufficiently close, validating the measurement [is validated] and 
calculating the average of the two curves [is calculated]. 
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20. (Amended) A calibration and nonnalization procedure, which are carried out 
[optionally] with the apparatus of claim [any one of the claims] 15[, 16,17 and 18 ], whereby the 
following steps are performed: 

1) in [In] order to perform relative intensity measurements, the differences of transmitted 
intensity between each fiber for the apparatus of claim 15 [or the differences of signal intensity 
between each optical detector for any of the claims 16 and 17,] are determined by performing a 
measurement on a turbid phantom illuminated uniformly or a diffusing sphere of perfectly 
uniform properties^. In this calibration procedure, the background light, measured with the light 
source turned off, must be subtracted from the signal. The obtained values are used to multiply 
the measurements given by each fiber and/or detector to correct the relative intensity 
measurements.] 

2) in [In] order to perform absolute intensity measurements, a caUbration [can be] performed on 
a turbid medium of known optical properties, which can be realized according to any one of the 
following recipes [recipe]: 

a) solid or liquid turbid medium which properties have been measured by other techniques, 
or reported in the literature, 

b) water suspension of micro-spheres of known size distribution and refractive index[. 
Absorbing dye may be added to the suspension. Li this case, the scattering properties are 
calculated from Mie theory, and the absorption coefficient is assumed to be equal to the water 
absorption coefficient, if no absorbing dye is added. If an absorbing dye is used, the absorption 
coefficient can be measured by a spectrophotometer, before mixing the solution with any 
scattering materials. 
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A Monte Carlo simulation is performed with the optical properties of the calibration sample. The 
simulation is then divided by the experimental reflectance performed on the caUbration sample. 
The result, that must be independent of the source-detector separation, is defined as the 
calibration factor. Each new measurement can be multiplied by the calibration factor]. 



-21 - 



17-11-200Q1 



CH 009900476 

RiC^dPCT/PTO 05 APR Ml 



09/806 831 



METHOD AND APPARATUS FOR MEASURING LOCALLY AND SIFPERFI- 1 

\ 

CIALLY THE SCATTERING AND ABSORPTION PROPERTIES OF TLTtBED 
MEDIA 



1 > Field of the invention 

The present invention relates to a method and an apparatus to quantify the optical 
scattering and absorption properties of turbid media, that can be applied in the extended 
optical domain of electromagnetic waves from far infrared (up to the microwave domain) 
to the far UV domain, i.e, where the energies of the photons composing the 
electromagnetic waves range from a few tnilli-eieclroa-Volt to more than ten electron- 
Volt More precisely the present invention relates to a non-invasive measurement, over a 
small area of the sample surface, Local and superficial characterization of biological 
tissues in vivo is a major application of this invendon. 

2. Related Background Art. 

Different techniques have already been proposed to quantitatively determine the 
absorption and reduced scattering coefficients of turbid medial Most of the non- invasive 
methods are based on the measurement of spatially and/or temporally-resolved 
reflectance- The principle is as follows: the turbid medium is illuminated by a coUimated 
or focused light source. The backscattered light is measured fay one or several detectors. 
Different types of measurements are possible, depending on the time-dependence of the 
illuminating source: steady-state (continuous source), rime-domain (short pulsed source) 

1 
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or frequency domain (amplitude modulated source). The present invention relates to the 
case of steady-state measurements, perfanaed at different distances /? between the source 
and the detectors. However, the technique presented here can be con^lcmented by tiroe- 
or frequBncy-domain measurernents. 

The range of values is an important point to consider, when comparing different 
methods based on the zneasurement of the reflectance. First, the probed volume of the 
turbid medium is related to the source-detector separation /?. The larger the source- 
detector s^aration, the deeper the average depth probed Second, depending on the range 
of A different mathematicaJ processing must be used to obtain the optical properties from 
the raw data. 

At least two cases must be distinguished 

1 ) The first case corresponds to source-detector separations f arger than several transport 
mean free paths. For typical biological tissue optical properties^ this case correspond to a 
source-detector separation larger than 2 mm. An analytical form of the reflectance can be 
obtained frora the diffusion equation, if the absorption coefficient fiQ is sufficiently lower 
than the scattenng coefficient /i, (typically ten times)- En such a case, the relevant optical 
propenies are the refractive index, the absorption coefficient and the reduced scattering 
coefficient. The average deptii of probing is on the same order than the source-detector 
separadon p. 



AMENDED SHEET 



1 7-1 1-2000 CH 009900476 



Such methods have been already published, and are the object of patents (Ref., patent 
5,517,987 Tsuchiya, Patent 5,676,142 Miwa et al,)- 

2) The second case corresponds to source-detector reparations close to one transport 
mean free path. For biological tissues', such source-detector separations correspond 
typically to distances fronx O.I to 2 mm. The average depih probed is on the order of 1 
mm. Such small source^detector separations enable the measurement of the optical 
properties of a small tissue volxnue. 

Wang et al. {Patent 5^630,423) proposed a method for the determination of the reduced 
scattering coefficient only, using an optical beam of oblique incidence. Moreover, their 
analysis does not include the effect of the phase function. Kessler et al. (Patents 
5,284,137 and 5,645,061) proposed a method for measuring the local dye concentration 
and scattering parameters in animal and human tissues, based on spatial and spectral 
measl^rements, Ho'^^ever, their methods do not enable the simultaneous deiermination of 
the absoTpt'ion coefficient, reduced scattering coefFicient. Other publications may be of 

concern; 

1 . See "Welch^ A. L; van Gemert, M. J. C. Optical Thermal Response of Laser Irradiated 
Tissue; Plenum publishing Corp., New York, 1995", and references therein. 

2. W,.F, Cheong, S.A. Piahl, and A J. Welch, "A Review of the Optical Properties of 
Bioiogica) Tissues/' IEEE L Quancuxn Electron. 26, 2166-2185 (1990). 
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3, Summary of the Inventioa: 

In the case of source-detecror separations close to one transporc mean free path, light 
propagation has been modeled using Monte Carlo simulations. In contrast to the previous 
case where the source-detector separation is larger than several transport mean free paths, 
it was discovered by the inventors that not only the first moment of the phase function 
must be considered, buc also the second moment More precisely, the relevant optical 
properties are the refractive index, the absorption, the reduced scattering coefficient and a 
new parameter I -g2)/(l-gi), that takes into account the first two moments of the phase 
function, denoted gi andgj respecrively. 

Different methods have been proposed for local characterization of turbid media, and in 
particular of biological tissues. Nfivenheless none allows for the simultaneous 
determination of the absorption, the reduced scattering coefSciem and the parameter 7. 

In the present invention we present a method and apparatus for the measurement of the 
absorption coefficient, the reduced scattering coefficient and the said phase function 
parameter 7, from the spatially-resolved reflectance data at short som'ce detector 
separation. These three parameters, which can be measured at different wavelengths, 
enable us to characterize turbid media, such as biological tissues. Our method is based on 
a theoretical study of the reflectance at short source-detector separations, that we 
performed with Monte Carlo siraulaiions. 
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BRIEF DESCRIPnON OF THE DRAWINGS 

Fig. 1 . Description of the principle of spatially-resolved reflectance measurement. 

/? denotes the source - detector distance, (Osoum and codtacwf are the solid angles of light 
acceptance determined by the Numerical Apertures (N.A-) of the optical fibers 

Fig.2 Probability densitj' function of the mean dqjth of scattering events per detected 
photons. 

Fig.3. Exampies of spatially resolved reflectance curves obtained with different phase 
funcrions. Case of matched refractive index Cn=I.O), 

Fig,4. Monte Carlo simulations and fits of the spatially resolved reflectance curves by 
mathematical expression in the formi?(/7) s j^/i (/?.//^,/) J . Case of constant 

and /i/- 1 am\ numerical aperture of the source and 
detectors = 0 37. A - 0.0647/-^^exp(.O.I6I/?), B ^ 0.18653 0,8466 mJ-^ 

Fig.S.a. Basic description of the apparatus; Case of measurements with an optical fiber 
probe. The apparatus comprises a probe coimected to a light source unit and a detection 
unit, both umt$ being connected to a computer or signal processing unit The optical fiber 
probe is constituted of one or more illuminating fibers, and one or more collecting fibers. 

S 
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Fig»5.b.Ba$ic descriptioii of the apparatus: Particular embodiment where the probe is 
composed of n deisctmg fibers (typically 6) m parallel. A broadband or white light source 
is used for excitation and the retro-di£fused light is coliected by &e n detecting fibers in 
parallel and dispersed in parallel in a spectrograph, before being detected on a 2D CCD 
camera and transmitted numericaily to a PC for jiMj'> Mc detennination at each 
wavelength^. 

Fig.S.c. Basic description of the apparaujs. Case of a measurement with sources and 
detectors directly in contact with the turbid medium. 

Fig-5.d. Basic description of the apparatus. Case of non-contact measurements with a ID 
or 2D detector, coupled with an imaging device. Scanning systems 1 and 2 are optional, 
they may be distinct from one another or conforaded in a single scaruxer. 

Fig.5-e. Basic description of the apparatus. Case of nOn-cOntact measurements with 
source and detection beams deflected by one or two scanning devices which may be 
distinct from one anoth^ or confounded. 

Fig. 6. Examples of the sample side of the fiber optical probe. 

6a. Simple arrangement for a single measurement 
6b. Arrangement for synunetrical measurements. 

6 
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6c. Arrangement for multiple measurements. 

6d. Airangement for multiple measurements, using an multicore fiber, or an optical fiber 
bundle. 

6e. Arrangement for symmetrical measurements where detecting fibers arc arranged on 
an elliptic path so as giving a regular spacing from the illuminating fibers to the detecting 
fibers, while maximizing the collecting surface of the fibers. Additional fibers F can be 
put in the center of the probe or around the probe for other light chzinnels: collecting 
fibers for fluorescence or Raman scattering for example. 

Fig. 7. Example of a measurement of the spatially-resolved reflectance with the probe 
shown in Fig.6.a. The measurement, performed on a microsphere suspension, is 
superimposed to a Monte Carlo simulation. The optical properties, computed from 
published water optical properties {jia) and Mie theory f^j' are; /i=1.33, ;i5=0.0004 
tmn\ ,u$-1^0 mm\ ?«2.2- The calibration was performed with a siloJtane sample of 
known optical properties. 



Fig.8, Relation between the parameters R(p-l mm) and \8 plnR(p^I mm}\ and the 
optical coefacients ftj' and pia .Case of /^Lf and 1.9. Probe of refractive index — 1.5^ 
sample of refractive index = 1.4> optical fibers diameter - 200 fino, NA-0.37 (source and 
collection). 
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Fig.9. Plot of the parameter ^>/^ for different /values (l.0. 1.8, 2.5), different reduced 

dp 

albedo a' (1, 0.95, 0.9) and for fixed p*^! mm. Mismatched reftactive index n«L4. 

Fig. 10. Plot of ^R(a'-l]-yjR{a') for different jrvalues and reduced albedo a', and for 
fixed p-l mm. Miamatched refractive index n«l .4. 
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1. Definitions 

The present definitions are given for clarity of the concepts developed in the patent and 
are worded in agreement with ref 1. 

The concept of spatially resolved reflectance is illustrated in Fig.h Consider a collimared 
light beam impinging on a turbid medium (through air or through a light guide, or even 
generated locally), in a given solid angle ©jource determined by the Numerical Aperture of 
the illuminating fiber: N.A. which is defined as the half angle of the entrance angle in the 
fiber, and on a limited area Ajt^urae^ where the source inadiate the surface. The spatially 
resolved reflectance R(p) is defined by the backscattered light power in a given solid 
angle cod^tecusr, at a distance p from the source, per unit area md normalized by the source 
power. The coHectmg area of the detector or detecting fiber is A44tfcicr. The distance pis 
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referred as the source-detector separsition (also in che esse when light guides or imaging 
devices are used for the illumination and for the collection of the backscattered H^t). 

R(p) depends on the optical properties of the turisid medium, defined below Note that 
R{p} depends also on die source and detectors chai3cieri5tics, i.e, numeri;;al aperture and 
siz^s. In a general situation, Ji{/:^ h the result of the convolution of die corrected 
reflectance curve and the source inadiance and detector sensitivity over, respectively, the 
surfaces A^oarce Md Ad6Kcior - When Asgurc^ and AdgtcciQ^ are smaller than the required spatial 
resolution > R(p) is confounded with the correaed spatially resolved reflectance. 

It is commonly admitted that the fundamental optical properties of a turbid mediunj 15 
determined by the average index of re&action n of the medium, the absorption coefficient 
Mi,, the scattering coefficient /^j, end the phase function p(0 where ^is the scattering 
angle. Tlie absorption coefficient /^a [rmn^] is defined as the probability of absarptioa per 
unit infinitesimal pathlength. The scattering coefficient /Zj [mm'^] is the scattering proba- 
bility per \mit infmitesinial pathlength. The phase function p( 0 is the density probability- 
function for the scattering angle The phase function is normalized as follow^: 



In turbid media and in biological tissues in particular, the phase function p^^is a 
continuous function of the scattering angle ^ and generally dqpends on the dielectric 
properties and on the material microstmcture. It can be approximated by a development 
in Legendre polynomials: 
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p{^hj-ti^"*^)sM<^o^^) (1.2) 



Where g« is the order moment of the phase function and P,^(^ is the Legendre 
polynomiai of order n: 

g.-2^[P„iff)p{ff)%inddff (1.3) 

The first moment of the phase function is also called the anisotropy fectar, and is often 
sinaply noted g (^^gj)^ It represents the mean cosine of the scattering angle. The reduced 
scattering coefficient jUs ' is defined as: 

ju/ ^ /u,fl-gj) (1.4) 
The transport mean fires path (or reduced mean free path) mfjp' of the photons, in a turbid 
medimn. is defined as: 

^^-f/W/W' (1.5) 

The reduced albedo a ' is the ratio: 

a'^Ms'/(M/^ M^) (1.6) 

As a result of the present invention described in more details in the next section, it is also 
necessary to define anotlier phase function parameter called 

r^fi'gMHs) (1-7) 
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All the parameters listed above are refcrxed as opUcaJ properties. They arc wavelength 
dependent, and can vary in space and time. In the following, /j gi and will be 
considered as parameters characterizing the nirbtd medium. 

2. Reflectance mea^tiremeats at distances close to one transport mean free p ath 

The naerbods described in die US Patent 5,517,987 (TsucMya) are based on 
measurements with large range of sonrce-detector separations, typically from I mSg' to 
10 mfp\ in suich cases, the volume probed is on the order of 10-1000 (mfp*)-*. In contrast 
with such a large scale investigation, the present invention is directed to a nove! approach 
where the volume probed is much smaller, on the order of 1 (m§)')^ This is achieved by 
using only small source detector separations, typically from 0.1 to 2 mfp\ The lateral 
dimension of probiag is limited to this range of distances. 

Xote that the present invention can be complemented by large source-detector separation 
measurements, if the optical properties of both superiiciai and deep parts of a turbid 
medium must be determined- 

A model of photon migration in tissue was necessary to predict the relationship between 
the measured reflectance and the optical properties. Analytical solutions from the 
diffusion equation are not appropriate in our case because we are interested in the 
reflectance close to the source, at distances comparable to the transport mean free path 
[mfp*]- We have performed Monte Carlo simulanons to predict the measured reflectance 
of an homogeneous semi-infinite turbid media. 

11 
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The exact diameter of the illuminating and collecting fibers, as well as flieir Numerical 
Apertures (N,A,X fae^^e been taken into account in the simulations. The mismatch of 
index of refraction at the surface of the medium have been also taken into account in the 
simulations^ by using the Fresnel law for each photon reaching the surface. 

This is also one result of this simulation to compute the average depth of probing, 
illustrated in Fig.2, . It is demonstrated thai only the superficial part of the turbid medium 
is probed if $mall source-detector separations arc used. For this, we determined the depth 
below the surface of each scattering event in the simulation. With this information we 
determined the average depth of all the scattering events for each detected photon, which 
we present as a probability density function in Fig.2. . This figure shows that the average 
depth of scattenng is approximately around I mfp'. Moreover it showed that the part 
located below 2 mfp' were not playing a significant role in the measured signal (for p< 
1.5 mfp-y 

The spatially resolved reflectance Rfj^, with short source-detector separations is' more 
complex than in the case of large source-detector separations. Indeed, for small source 
detector-separation conditions, the mverse problem, i.e. calculating the localized 
absorption and reduced scattering coefficients, is sensitive to the scattering phase 
mnction. It is part of the present invention to have shown, from Monte Carlo simulations, 
that only the f5rst and second moments of the phase function can be taken into account. 
Moreover, it was established that the influence of these two moments are not 
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indepejodent. Indeed, it is a merit of this invention to demonstrate that only one parameter 
p^C^'g^^fJ'giJ* which depends on the first and second moments of the phase fimction 
can be used to characterize acciirately the reflectance at short source-detector separation. 

Fig.3 illustrates the fact chat the parameter /is the only predominant parameter of the 
phase function that must be taken into account (and not g^gj as frequently mentioned in 
literature). Different leflectance curves, obtamed from Monte Carlo singulations are 
shown in Fig.3, Four different phase functions \vere lised for the simulations. Three phase 
functions are characterized by /-1 25, but different and values, Fig.3. shows that 
aLcnost identical reflectance curve are obtained for distances /yij ' >0.3, if the gj and g-^are 
varied in such a way that the parameter 7^ stays constant (;^-L25). For comparison, Uie 
reflectance computed wirh a different g{'='§0 value is also presented (r=2.25). It appears 
then clearly that the paranoeter g{=gi) plays a significant role in the reflectance. 

The pararaeter /depends on the characteristics of the set of scatterers inside the turbid 
medium: shape^ distribution of sizes and refractive index. 

The possibility of determining this parameter ;^is an important part of this invention, and 
has never been reported earlier. The determination of ^f'can be related to the microscopic 
strjcture of the sample, since y is related to the sizes and refractive mdexes of the 
scatterers. 
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In most cases, tbe average refractive indices n and die soiirce and detectors characteristics 
are fixed and known. In such case, only three parameters are determined from measure- 
ments at short source-detector separation: iia, ' and ;< 

It is important to note that the parameter cannot be estimated if large source detector 
separations are iised. Indeed, the reflectance is only sensitive to ^^at short source-detector 



It was also derived from Monte Carlo simulations that tte reflectance R(p), for Imfjp', 
can be reasonably well approximated by the expression: 



it is important for the description of the invention to evidence the fact that the fimction A 
depends on a the scattering properties (i.e. ' )) but not on ^a- ^ contrast, tiie ftaxc- 
tion B depends on and ii^ ' but neither on >'Tior p. An example of Equ. (L9) is shown in 
Fig. 4. 

This formulation is of great help to solve the inverse problem, as described in section 5.3. 



separation. 




(1.8) 



In many cases, we foTuid that Equ.( 1 . 8) can also be written as: 




(1.9) 
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3. Apparatus 

In the first embodiment. The apparacus can be divided in three parts, described in Fig.5.a. 

The first part i$ the illuxniMting system. Any light source can be used. For examples: a) 
white sciiTces, such as halogen or xenon lights, metal halides or fluorescent or 
phosphorescent sources, b) sources such as lasers, laser diodes, optical fiber lasers, light 
emitting diodes or superluminescent diodes, c) the sources described in point a) and b) 
where monochromators, filters or interference filters are added to select a given set of 
wavelengths. 

In the first preferred embodiment, The light power is conducted to die investigated 
sample by the probe, which is d:e second part of the apparatus. The probe is preferably 
made of optical fibers, to illuminate and to collect the backscattered light But GRIN rods 
or other type of light pipes can also be used. Different possible anangements of optical 
fibers are illustrated in Fig-6 . Two different modes of measurements can be chosen. First, 
one fiber is used to iiluniinate the sample and at least two other fibers are used to collect 
the backscattered light at two different distances. Second, one fiber is used to collect the 
light and at Itast two other fibers, located at two different distances from the first one, are 
used to illuminate sequentially the sample. 

The arrangement of the different flbei^ can be replaced by any imaging system or image 
guide, such as multicore optical fibers, or optical fiber bundle. 
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The light collected by the probe is annlyzed by the detection unit, which :s the third part 
of th& apparatus. If wide spectral light sources are used (such as halosen or xenon tights), 
a spectrograph can be put between the probe and the detector to get wavelength 
dependence of the backscattered signal (either in the source or detection unit). Any types 
of detectors can be used. For example, photodiodes, avalanche photodiodes or 
photomultipliers can be assigned to each collecting fibers. Simultaneous detection of each 
collecting fiber can also be achieved using Linear or two-dimensional detectors such as 
Charge-Coupled Detectors (ID or 2D), intensified CCD or array of photodiodes, 

A particular embodiment is described m Fig.S.b. : the probe is composed of n detecting 
fibers (typically 6). A broadband or white light source is used for excitation and the retro-* 
diffused light is collected by the n detecting fibers and dispersed in parallel in a 
spectrograph, before being detected by a 2D CCD camera and transmitted numerically to 
a PC for \ fia determination at each wavelength X, This particular embodiment has 
the advantage of yielding the wavelength dependence of from a single 

measurement. 

A second type of embodiment is presented on Fig. 5.c. The difference with embodiments 
presented in Fig. 5. a. and Fig. 5.b. is that optionally no optical fibers, light pipes or grin 
rods is used The light source unit is directly in contact with the turbid medium, as well as 
the detector imit. Coilimating optics, micro-optics or imaging optics (DOE Dif&active 
Optical Elements for example) can be put between the turbid medium and the actual light 
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sources and detectors. The different type of sources and detectors cited in example for the 
first embodiment can be used for ihe second type of embodiments. Hybrid design, such as 
airangements involving both direct contacts sensors or detectors and fibers, light pipes or 
grin rods are also included in the present embodiment. 

The third emfaodimeni is descnbed in Fig.5.d and Fig.S.e. The only difference with the 
first and second embodiment is that non-contact measurements are performed, A 
coUimating system allows for a point-like illumination on the turbid medium. An imaging 
system enables the measurement of the spatial distribution of the reflectance. The 
detectors can be either an array (ID or 2D) of detectors (Fig.S.d) with optional scanning 
systems 1 and 2 which can operate separately or can be confounded in a single scanner, 
or a single detector (Fig.S.e). In this last case, a scanning device is used to obtain the 
spatially-resolved reflectance. A fiber bundle, multicore fiber or relay optics (grin rod or 
multiple lenses) can be put between the focal point of the imaging system and the 
detecior($). The different type of sources and detectors cited in example for the first 
embodiment can be used for the third embodiment 

Spatial images of the parameters >u/ and / can be obtained by a series of 
measurements at different locations; that we call multi»site measurements. Some 
mechanical or optical scanning device can be used for this purpose. The resolution of 
such images is on the order of the mean soxirce-decector separation used for a single site 
measurement. All three embodiments can be e?cpanded to perform multi-site 
measurements, by duplicating and/or multiplexing the illuminating or measuring devices. 
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For example, ±e optical fiber probes shown in Fig.e^ can be duplicated and put side by 
side. Tile scanning system described in the third embodiment cai also be expanded to 
perform multi-site measurements. 

4. Normalization and calibration 

The differences of transmission between each fiber are corrected by performing a mea- 
surement on a turbid phantom illuminated uniformly. 

The background li^t, measured with the light source turned off, must be subtracted from 
the signal 

In order to perfonn absolute intensity Eneasurements, a calibration is performed on a 
turbid medmm of known optical propercies. Examples of such a medium are: 1) solid or 
liquid turbid medium which properties have been measured by other techniques, or 
reported in the literature, 2) water suspension of microsphere of known size distribution 
and refractive index. Absorbing dye may be added to the suspension. In case 2) the 
scattering properties are calculated from Mie theory, and the absorption coefficient is 
assimed to be equal to the water absorption coefficient, if no absorbing dye is added. If 
an absorbing dye is used, the absorption coefficient can be measured by a 
spectrophotometer, before mixing the solution with any scattering materials. 

A Monte Carlo simulation is performed with the optical properties of the cah1>ration 
sample. The simulation is then divided by the experimental reflectance performed on the 
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calibration sample. The result, that must be independent of the source-detector separation, 
is defined as the calibration fector Each new measwement is multiplied by the 
calibiation factor. 

An example of a measxirement, after calibration, i$ shown in Fig J. The measurement was 
performed with the apparatus described in the first embodiment, with a probe similar to 
the one described in Fig.6,a. A prior calibration was performed on a solid turi^id medium, 
which optical propcnies were measured by another technique (frequency domain 
measurement, cf, Ref L). The sample was tt water suspension of pol>^tyrene microsphere 
of I Jim diameter. The measurement is compared to a Monte Carlo simulation performed 
with the scattering properties calculated from Mie theory and absorption coefficient of 
water, shows an excellent agreement between the measurement and tihe simulation, 
which confirm ths validity of our simuladon model, experimental measurement and 
calibration. 

5._Signal processing 

5.1. Contrpl of the homogeneity of the area probed 

Artifacts during a measurement, for example due to bad contact between The probe and 
the sample, or heterogeneity of the sarople, can be detected by the following optional 
procedure. Two illuminating fibers are disposed s)TOmetricaUy in regard to the collecting 
fibers (see Fig,6.b). If the sample is homogeneous, the reflectance curve should be 
identical with either illuminating fiber. Therefore, heterogeneity of the investigated 
region or obstructions beneath the fibers are detected by comparing the two curves. If the 
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two curves are sufficiently close, the measurement is validated and the average of the two 
curves i$ calculated. 

5 .2 . Saioothmg procedare and computation of the derivative of &e curve . 

Functions in the form m^p"^^ ^p[m^p) were always found to fit well Monte Carlo 
smiulatioGS for a restricted range of distances. Therefore, a smoothing of the 
experittiental reflectance R(p) con be proces^d by fitting R(pJ to m^p"^ ^vij^-^P) - 

The determination of the slope of the logarithm — (in/t(/?,/jj,^',;^)) is also derived 
from this fit, using the following fonnuia: 




(IJO) 



The slope of the square root of R(p} can also be derived (the use of these slopes are 
described below in section 53 ): 




(1 H) 



53. Inverse problem 

The so-caJled **direct problem*' consists in computing the spatially resolved reflectance 
Kip) from the values of the optical coefEcieats /^j fjis\ y involved in a model of 
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propagation of die tight in turbid medium and whereby the said "model" incorporates a 
Legendre polynomial deveiopment to the second OTder of the said '*phase function", and 
whereby the said "phase parameter" y is introduced in the computation as an 
independent parameter. The so-called "inverse problem" consisis in extracting die optical 
coefficients /Xfl,.^j* 7" from the spatially resolved reflectance data . 

We developed different methods to solve the inveise problem, which consists tn 
extracting optical coefficients ftom the reflectance data, 

Method 1. 

Monte Carlo simulations show that the measurements of the reflectance intensity R(p) 
and the slope of InAf"^ (denoted dplnR), determined at a fixed distance /?, can be used 
to derive fij ' and for a given ;>'value. Fig.8> shows graphically the relationship between 
and jia and the two parameters R(p-1 mm) and \ dplnR(p^l mm)\ in the case of 
p^l mm, L5 and 1.9. Note that any other choice of p\% possible, as long as ^ is on 
the order of 1 mfp'. We see in Fig.8. that ' and jia can not be determined uniquely ify 
is unknown, further insight for optimising the inversion strategy is provided by diree 
additional features of 

First, the determination of ' is only weakly influenced by y Indeed in Fzg.8. the errors 
induced by error in /are typically 10% for/i/ around 1 mm"^ Second, although absolute 
detenninadoG of is not possible when y is not precisely known, relative absorption 
changes can be still evaluated. Third, the indetermination of the parameter / may be 
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resolved by the values of R(pJ and/or | dfinR(p)\ at other distances. Therefore the follow, 
ing procedure can be used: 

{ 1) determination oifXs ' and Ma from R{p ^1 mm) and \dplnR(p «i mm)\ for a set of 
values ;^ For example: ;k=1.0, LI, 1.2,...,2.5. 

(2) simulations with the different sets of fi^ ' and }Xa obtained from point 1) 

(3) comparison between the simulations and the reflectance profile for distances 
0.3</><2 mm. 

This last step allows us to detennine the correct values ofy, fis' and ^la- Points 1 to 3 can 
be done iteratively to evaluate rinore precisely, using a finer discrimination of /values. 

This method described by the iteration of points 1) to 3) corresponds to finding a 
simulation that fit closely a measured curve Rfp). This can also be achieved by directly 
comparing the measured curve Rfp), with a set of simulations, and finding the simulation 
that minimize the square of the differences between the measured and simulated curves. 

Note also that the use of the reflectance intensity R(p) and the slope of !nR{p) is 
equivalent to the use of two intensities R(pi) and Rfpzl. measured at two close distances 
Pi and p2, respectively. 
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Method 2r 

The inverse problem can also be performed^ considering properties of i?^^ expressed by 
Equ.(L8). 

Indeed, it can be derived from Equ.(1.9) diar flie quantity docs not depend on the 
absorption coefficient fia: 



j^^R{p,M„K,r) ^^{p,f^,,r) (1.12) 



g I ■ ■ ■ 

This property i$ confrnned in Fig.9 , where the quantity —AR{p^^^,p^^y\ , evaluated 

at /7"1 mm, is plotted as a function of for , 1 .9 and 2.5 and reduced albedo a -1, 
0-95 and 0.9 (note that any other choice of p is possible, as long as p is of the order of I 
mfp^). 



In Fig.9. , the parameter —jR{p,p:^,f:^.y) clearly depends onf^/ and In contrast, 
the dependence on a' is almost negligible. 

Therefore, /and ji/ can be derived from the parameter —^Ri^p^.fr^^p^,}''^, calculated 

from the experimental reflectance R(p) (see section 5.2). Simultaneous determination of y 

and Us ' require values of —JR(pjP'^,p^,y)ax different distances p (at least two). If one 

dp ' ^ ' 

of the parameters /or /x,' is already known, the detennination of the other one can be 
obtained S-om the value of -^-^^(A/^./^a./) ^ single distance. 
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Convenient analytic^ approximation of — Jiff Ai^CtAi,./) can be obtained by fitting 

op ^ 

Monte Carlo results to polynomial functions- For example, in the case of fixed /^1*9, we 
obcained: 

Ms '-0.9162-S2,8»x^l795x^.l8l55x^+65428x* (lU) 



where x = --^^R{p,/r,.A./) at p=l mm. x is expressed in [mm*^] and is 
expressed In [mm'^]. 

This results is valid for a probe of refi^tive index of 1.5 and a sampte of refractive index 
of 1.4, optical fibers of diameter 200 |im, NA=0.37 (source and collection). 

Once Ms ' and /are calculated from the procedure explained above, /ig is calculated j&om 
the absolute value of R(p), which highly depends on p^. F or given jUs \ ;'aiid p values, the 
dependence of reflectance on fia is obtained by Monte Carlo results. From Equ.(1.8) we 
have: 

where and are functions given by Monte Carlo simulation. Panicularly, they can be well 
approximated by polynomial functions. For example, for ;*1.9, probe refractive index of 
1.5, sample refractive index of 1.4, optical fibers of diameter 200 \im, NA^J7 (source 
and collection): 

h - -0.002257 - 8,171 y + 268.8 y (1.15) 
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/= 0,0131 RO,051S4 /As -0.01974 pt/ "V 

0,003217 /£/ ^-0.0001992 /z/ ^ (1.16) 



where / = 



^R(p,^,.Mc*y) -/][mm'^]. //sis expressed in [mm'']. 



M:etho_d3 

Equ.{l,8) also shows that relative measurements of fia^ \x, variation of /^a from a known 
value ^iffo, is possible by measuring the variation of the parameter ^ J! ( a ^ * A, , r ) ' 



This reladon is illustrated in Fig.lO. in the ca$e of /x^^^O, /?=J mm, ;r==LO, L9 and 2.5, 
and for a' == 1 to 0.83, confinns that the influence of ;^ is weak in the 



quantity - ^R^p.M'.^J^a^f) - known p./, the function B(a') allows 
for. the determination of a relative absorption change dfia "/iff - ^<n>» Fig. 10. illustrates the 
case Mao -0, but any other value of /u^o is possible. The interesting point is that B(a ') does 
not depend on the phase function. 
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CLAIMS 

A method for local and superficiaJ (probed volume and meastiiing area with 

dimensions on the order of one transport mean free path) characterization of a 

turbid medium using the following parameters: 

1} the refractive index n of the turbid medium 

2) tbe absorption coefficient jx^ of tiie turbid medium, 

3} the reduced scattering coefficient /i^' of the turbid medium. 

4) the phase function parameter r=0'S^^O*S'') of the turbid media, where gj and 

g2 are the first two moments of the Legendre polynomial development of the 

phase function p, 
and comprising the steps of: 

-measunng the spatially-resolved reflectance Rfpjofihc turbid medium 0? being 
the source-detector distance) by any mean^ comprising an {Uumination beam as a 
source and an optical detector, which, by using optional signal processing, which 
may involve filtering and de-convolution operation to correct for the non-zero 
area of either the illumination source or the detector, allows for the precise 
determination of the said spatially-resolved reflectance R(pJ, 

- marhematically processing Rfpj to compute at least one of die said parameters: n, 
r,Ms\ r and/or the variations, in time and/or space, of at least one of the said 
parameters: An, A/i^ A^^^', Ay , whereby the said "inverse problem", which 
consists in extracting the optical coefficients from the spatially resolved 
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reflectance data is solved, and whereby the said "direct problem" cojisists in 
computing tiie spatially resolved reflectance from the values of the optical 
coefficients n. Ma Ms '> r involved in a model of propagation of the light in turbid 
medium and whereby the said "model" incorporates a Legendre polynomial 
development to ±e second order of the said "phase function" and whereby the 
said "phase parameter" r is introduced in the computation as an independent 
parameter. 

The method of claim I , wherein said spatially resolved reflectance is measured by a 
probe comprising at least one optical fiber canying the light from the source unit to 
die mrbid medium and at least one optical fiber collecting the reflected light and 
canying it to detection unit, whereby the combination of a variety of emitting ftbsrs 
and of receiving fibers yields a set of distances p at which the reflectance R{p) is 
measured. 

The method according to any one of claims 1 and 2, wherdn said spatially resolved 
reflectance R(p) is measured for a set of values ofp, by using a probe composed of 
optical fibers in any of the following configurations; 

- one emitting optical fiber and a set of optical receiving fibers 

- a set of optical emitting fibers and one optical receiving fiber 

- a set of optical emitting fibers and a set of optical receiving fibers 
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giving the spatially resolved reflectance R(f^ at a varieiy of source - detector 
distance p. The emitting and receiving fibers can be arranged in a variety of 
configurations. In some particular arrangement they can be positioned: 

- on a line, 

- on crossed Une$» 
• on a circle, 

- on an ellipse 

- on crossed ellipses 

" on a diskj a rectangle, or any surface, as a dense and contiguous arrangement of 
fibere, such as it can be obtained with a bundle of fibers or a tnulti-core fiber, 

- on any pattern resulting from the combination of the above mentioned patterns. 

The method according to any one of claims 1, 2 and 3, wherem said spatially 
resolved reflectance is measured by an oprical and electronic micto-system 
comprising a coUimated or focused beam as illuminating source and ID or 2D 
arrays of optical detectors, arranged in a way similar to the one described in claim 
3, and which, in some embodiment, can be a CCD or MOS camera. 

The method according to any ooe of claims 1-4, where the probe of claims 2 and 3 
or the optical and electronic micro-system of claim 4 are put in contact to the turbid 
medium. 
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6. The method of claim 1 , wherein said spatially resolved reflectancs is measured by a 
non-contact system, comprising at least one of the following combination of optical 
systems: 

- a fixed optical system to irradiate the turbid medium with a coilimated or focused 
beam forming the illuminating source and a fixed optical system comprismg an 
imaging system forming the image of the measured area of the turbid medium on 
a said "optical detector^, which can be foreaed of ID or 2D array of optica! 
detectors, whereby this second optical system can be identical to the first one and 
whereby the array of optical detectors can be either one of the following system: 

• a set of optical fibers arranged in a variety of configurations similar to 

those described in claim 3, 
■ an optical and electronic micro-system (MOEM), which in some 

embodiment can be a CCD or MOS camera, with a functionality similar to 

the one described in claim 4, 

- a fixed optical system for the coilimated beam illuminating source and a scanned 
optical system for the said "optical detector*', 

- a scanning optical system for the coilimated beam illuminatiag source and a fixed 
optical syscem for the said ''optical detector**, 

- a scanning optical system for the coilimated or focused beam used as an 
illuminatins source and a scanning optical system for the said "optical detector**. 

7- The method according to any one of claims 1. to 6., wherein the absoiption 
coeiScient ^a, the reduced scattering coefficient jis' and the phase function 
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parameter 7 are detennined by fitdng the me^ured spatially-resaived reflectance 
^(fl M/f Mi» Ji^ to « discretized data obtained by using Monte Carlo 

simulations, or to interpolating functions giving a continuous approximation of the 
discretized data obtained by Monte Carlo simulations, and whereby said "Monte 
Carlo simulations'' are based on a photon propagation model comprising a phase 
ftinciion approximated by a Legendre polynomial development limited to the 
second order. In a preferred embodiment, cubic splines are u$ed as interpolating 
functions. 

8. The method according to any one of claims 1 to 7, wherem the following signal 
pToce$$mg steps can be performed in addition to the processing steps of any of the 
claims 1-6 whereby the processing of the spatially resolved reflectance data of 
claim 7 can be simplified and accelerated by: 
- fitting the measured reflectance R{p} by the function: 

This fit gives the values of the parsmeters mi, and 013 , assuming that the 
expression R(p)=m^p'"' exp(^m^p) gives a smoothed expression of the spatially 

resolved reflectance Rfp). 

-Computing the slopes —J Rip) and ~{\nR{p))^ or any mathematical 

dp dp^ ^ ^' 

combinations of these two latter quantities and R(p), from analytical functions 
using the parameters mi. nx-^, or by numerical proceduuss from the expression 

R(p)^m,p'^'ex:p{m^p) 
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- computing the values of at least one of Ae said parameters: «. ,. f^s r snd/or 
the variations, in tiroe and/or space, of at least one of the said parameters: An, 
^tfg, Ciji:,'. iky from the relationship between the reflectance intensity R(fi^ and 
the slope of In/fC/jl (denoted dplnR), detennined at a fixed distance 
/? comparable to the transport mean ftee path, whereby the computation is made 
from the data obtained by Monte Carlo simulations, provided that said "Monte 
Carlo simulations" are based on a photon propagation model comprising a phase 
function approximated by a Legendre polynomial development limited to the 
second order. 

The method according to any one of claims 1 to 7, wherein the following signal 
processing steps can be perfoiiiied in addition to the processing steps of any claim 
1-6 whereby the processing of the spatially r^olved reflectance data of claim 7 can 
be simplified and accelerated by: 

- computing the reduced scattering coefficient ^/ and the phase function parameter 

by using the following form of Hxe reflectance: 

RfpJ = (ACa r. (is ') fi, 

where the function A(/;i % t^^') and Bfpta, fi,') depend also on the sources and 
detectors chart.cteristics, and the refractive index of the sample, and comprising 
the steps of: 
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- Computing the slopes of ihc square root of the $patially resolved reflectance 



coefficient for Q3<f^s '<5, for at least two distances />, 
- Determining the parameters and / by a polynomial mteipolarion of the data 
obtained by Mome Carlo simulations, whereby $aid "Monte Carlo simul^ons" 
are based on a photon propagadon model comprising a phase function 
approximated by a L^endre polynomial development limited to the second order. 

10. The method of claims 9, wherein the absoiption coefficient pta is determined by 
using the equation: 



where /and h are continuous functions of the parameters ju/ and ^of that can 
be approximated by a polynomial interpolation of the data obtained by Monte Carlo 
simulations, whereby said '*Monte Carlo simulations" are based on a photon 
propagation model comprising a phase function ^proximated by a Legendre 
polynomial development limited to the second order, 

11* The method according to any one of claims 1 to 7, wherein the difference 
^Ua-^a- Maa between the absorption coefficient p.^ and a known vajue is 

computed from the quantity ^i?(/?,^,/^^o.r)->/i2(A/<;»Afl»r) = B(a')' B(ao% 
whereby tiie function B(a') of the albedo a' can be can be approximated by a 



=r-(/3,/<;.r), which d^ends weakly on Hie absorption 
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polynomial inteipolation of the data obtained by Moute Carlo simulations, whereby 
said ''Monte Carlo simulations" are based on a photon propagation model 
comprising a phase function approximated by a Legendre poiynomiai development 
limited to the second order^ and whereby the calculation can be done for a single / 
value, because the influence of the phase function parameter and j^in B(a') are 
pardculariy weak. 

The method according lo any one of claims 1 to 1 1, wherein the illuminating source 
is a broadband source or a white light source and the detector unit comprises a 
spectrograph or any wavelength analysis system to measure the wavelength 
dependence of at least one of the parameters (n, pa^jij \ . 

The method according to any one of claims 1 to 12, wherein said turbid medium is 
a biological medium. 

The method according to any one of claims 1 to 13, wherein the measurement and 
processing is repeated at different location of the sample (multi-site measurements), 
to build images of any one of the said parameters (n, fix^s\ 

15. An apparatus using the method of claim i to 3, 5 or 7 to 1 4 for local ajxd superficial 
characterization of a turbid medium, 

a) comprising a source, a detection unit> reference means, signal processing means, 
a probe comprising at least one optical fiber connecting said source unit to the 
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turbid medium and at least one optical fiber connecting the turbid medium to the 
said detection urdtj and ref«"ence means 
b) where the distance between the source and the detector is close to one transport 
mean free path 

16, An apparatus using the method of claim I, 4, 5, or 7 to 14 for local and superficial 
characterization of a turbid medium, 

a) comprising an optical and electronic micro-system comprising at least one 
illuminating source, at least one detector , signal processing means and reference 
means ^ 

b) where the distance between the so mice and the detector is close to one transport 
mean &ee path 

17, An apparatus using the method of claim 1 or 6 to 14 for local and superficial 
characterisation of a mrbid medtimi, 

a) comprising a coUimated or focused beam used as an illuminating sotirce, at least 
an optical detector for the detection unit, a fixed or scanning optical system for the 
illimiinating source and a fixed or scanning optical system for the said "optical 
detector", signal processing means and reference means, 

b) where die distance between the source and the detector is close to one transport 
mean free path 
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18, Three apparatus similar to the apparatus described in any one of the claims 15, 16 
and 17 , and characterized by the feet that the distance between the colltmated or 
focused optical beam used as iiluminadng source and the emitting point connected 
to an optical detector varies &om 0.1 to Zmm. for application to biologicaJ media 
and to turbid media having a transport mean 5:ee path similar to biological media. 

19, A test , where the control of the homogeneity of the sample over the probed area is 
performed v^^ith the apparatus of any one of fiie claims 15, 16,1? and 1 8 ♦ which can 
be carried out according to the following optional procedure: In die apparams two 
illuminating fibers aie disposed symmetrically in regard to the collecting fibers. If 
the sample is homogeneous, the reflectance curve should be idenrical with cither 
illuminating fiber. Therefore, heterogeneity of the investigated region or 
obstructions beneath the fibers are detected by comparing the two curves. If the two 
curves are sufficiently close, the measurement is validated and the average of the 
two curves is calculated. 

20, A cah'bration and uormalization procedure, which are carried out optionally with 
the apparatus of any one of the claims 15, 16,17 and 18 , whereby the following 
steps are performed: 

1) In order to perfoim relative intensity measurements, the differences of 
tninsmitted intensity between each fiber for the apparatus of claim 1 5 cr the 
differences of signal intensity between each optical detector for any of the 
claims 16 and 17, are determined by performing a measurement on a turbid 
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phantom Ulumniated uniformly or a difiusiisg sphere of perfectly uniform 
properties. In this calibrstion procedure, the background light, measured with the 
light source turned off, must be subtracted from the signal. The obtained values 
are used to multiply the measurements given by each fiber and/or detector to 
correct the relative intensity measurements. 
2) La order to perform absolute intensity measurements, a calibration can be 
f^ormed on a turbid medium of known optical properties, which can be 
realized according to any one of the following recipe: 

a) solid or liquid turbid medium which properties have been measured fay other 
techniques, or reported m the literature, 

b) ^vater suspension of micro-spheres of known size distribution and refractive 
index. Absorbing dye may be added to the suspension. In this case, the 
scanering properties are calculated from Mie theory, and the absorption 
coefficient is assumed to be equal to the water absorption coefficient^ if no 
absorbing dye is added. If an absorbing dye is used, the absorption coefficient 
can be measured by a spectrophotometer, before mixing the solution with any 
scattering materials. 

A Monte Carlo simulation is performed with the optical properties of the 
calibration sample. The simulation is then divided by the experimental reflectance 
performed on the calibrarion sample. The result that must be independent of the 
source-detector separation, is defined as the calibration factor. Each new 
measurement can be multiplied by the calibration factor. 
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